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Abstract 


Solar activities have a great impact on modern high-tech systems, such as human aerospace activities, satellite 
communication and navigation, deep space exploration, and related scientific research. Therefore, studying the 
long-term evolution trend of solar activity and accurately predicting the future solar cycles are highly anticipated. 
Based on the wavelet transform and empirical function fitting of the longest recorded data of the annual average 
relative sunspot number (ASN) series of 323 yr to date, this work decisively verifies the existence of the solar 
century cycles and confirms that its length is about 104.0 yr, and the magnitude has a slightly increasing trend on 
the timescale of several hundred years. Based on this long-term evolutionary trend, we predict solar cycles 25 and 
26 by using phase similar prediction methods. As for solar cycle 25, its maximum ASN will be about 
146.7 + 33.40, obviously stronger than solar cycle 24. The peak year will occur approximately in 2024, and the 
period will be about 11 + 1 yr. As for solar cycle 26, it will start around 2030, and reach its maximum between 
2035 and 2036, with maximum ASN of about 133.0 + 3.200, and period of about 10 yr. 
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1. Introduction 


The Sun, as the unique star in the solar system that is closest to 
the Earth, plays a crucial role in the formation and evolution of the 
Earth and its near-Earth space environment. At the same time, 
solar activities may produce great damage to modern high-tech 
human systems, such as aerospace, satellite communication and 
navigation, large power grids, deep-space exploration, national 
safety and related scientific research. Studying the long-term 
trends in solar activity will help us answer the big questions such 
as how the Earth’s future environment will change and where 
humans will go. At least, it can also help us predict the intensity 
and development trend of future solar activity. 

As is well known, solar activity has cyclical characteristics, 
reflected in various indicators, such as sunspot number, sunspot 
area, 10.7 cm—wavelength solar radio flux, etc. One of the 
most obvious periodicities is the solar cycle with an average 
period of about 11 yr, also known as the Schwabe cycle, first 
reported by Schwabe & Schwabe Herrn (1844). 

Over the past 200 yr since the discovery of the Schwabe 
solar cycle, some evidence with other periodic characteristics 
has also been discovered, such as the grand minimum, the G-O 
rules that show the difference between even and odd cycles and 
the century cycle (Hathaway 2015). At present, some of these 
features are well understood, but others are not yet conclusive. 
For example, does the century cycle really exist? How long is 
its period? Does it change? If so, how does it change? 


The solar century cycle refers to a periodicity of approximately 
around one century, first proposed by Gleissberg (1939). He 
smoothed the amplitude of sunspot numbers occurring between 
1750 and 1928 and found a long period of 7-8 solar cycles, which 
is close to a century and also known as the secular cycle or 
Gleissberg cycle. In 1994, Rozelot (1994) used Fourier analysis 
and proposed that the period of the century cycle is 97.2 yr. 
However, Garcia and Mouradian’s analysis (1998) four years later 
suggested that the period of the century cycle is 78 or 81 yr. Other 
researchers (Peristykh & Damon 2003) obtained evidence that the 
period should be 88 yr by using the analysis of cosmogenic 
isotopes. Ogurtsov et al. (2002) more explicitly described the 
century cycle as a mixture of biperiodic structures with a 50-80 yr 
period and a 90-140 yr period. Le & Wang (2003) proposed that 
the period should be about 101 yr. Recently, some researchers 
(Ma 2009) have extended the length of the century cycle to 
between 60 and 150 yr. Tan (2011) analyzed average relative 
sunspot number (ASN) during 1700 AD-2009 AD and proposed 
that there should be three kinds of solar cycles: the well-known 
11 yr Schwabe cycle, the 103 yr century cycle and the 51.5 yr 
cycle. The periodicity of the century cycle also appears in solar 
proton events and is possibly linked to the grand minimum 
(McCracken et al. 2001). 

The causes of the century cycle were analyzed. In 1999, 
Pipin (1999) proposed a numerical model of the century cycle 
based on the dynamo mechanism describing the generating 
processes of solar magnetic fields, and estimated that the 
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century cycle is the result of the re-establishment of differential 
rotation after the magnetic feedback of angular momentum 
transport. Some researchers take a different view of the century 
cycle. In 1999, Hathaway et al. (1999) used the century cycle to 
predict future cycles and found that the best-fitting value of the 
century cycle was constantly changing rather than stable. 

Schwabe (1843) made the first solar cycle prediction (Attia 
et al. 2020). After almost two hundred years of development, 
Petrovay (2010) divided solar cycle prediction methods into 
three categories: methods based on physical models, experi- 
ence-based precursors and extrapolation methods. 

The model-based method mainly refers to the dynamo model of 
the solar magnetic field. The precursor method relies on previous 
solar activity indicators or magnetic fields at earlier times to 
predict the maximum amplitude of the next solar cycle 
(Petrovay 2020). It is mainly divided into two categories: 
geomagnetic index prediction and polar magnetic field prediction. 
The prediction of the polar magnetic field as a precursor was first 
proposed by Schatten et al. (1978). Nandy (2021) believes that the 
precursor method is relatively superior to other methods in terms 
of accuracy, especially the approach that uses the polar field near 
the solar cycle minimum as the precursor factor. The use of the 
geomagnetic index as a precursor dates back to 1966. In that year, 
Ohl found that the value of the aa index indicating geomagnetic 
activity near the time of the solar cycle minimum correlated well 
with the amplitude of the next cycle (Hathaway 2015). According 
to the definition of the precursor method, not only can the 
geomagnetic index and the magnitude of the polar magnetic field 
be utilized as precursors, but the sunspot numbers in the declining 
phase of the cycle (Brajša et al. 2022; Nagovitsyn & Ivanov 2023), 
the number of spotless days in the declining phase of the cycle 
(Burud et al. 2021), etc. can also be regarded as precursors. 

The extrapolation methods are based on statistical analysis of 
various indicators of solar activity to model its behavior, then this 
model is applied to reasonably extrapolate and predict the 
following solar cycles. Recently, this method was used to predict 
solar cycle 25 (Kakad & Kakad 2021). There are many metrics 
that the extrapolation method can utilize, but the most widely used 
is the sunspot number sequence. Sunspots have been observed for 
more than three centuries and can be regarded as the longest 
running time series in the world (Bhowmik & Nandy 2018), 
which is one of the main reasons why the sunspot sequence was 
chosen as a predictor in this work. Additionally, similarity theory 
also can be used to predict the forthcoming cycles which assumes 
that two solar cycles with similar solar activity indicators should 
also have similar features (Du 2023). 

In addition to the above classification, cross-disciplinary 
approaches such as machine learning and neural networks for 
solar cycle prediction are also very popular research directions. 
Ramos et al. (2023) summarize the results on using machine 
learning to make predictions for solar cycle 25. 

The methodology of this work is based on the extrapolation 
method and similarity theory. First, we confirm the existence 
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and long-term evolution trend of the solar century cycle, 
obtaining an empirical century cycle function by using the 
fitting method. Then we propose that an 11 yr solar cycle at a 
similar phase in the century cycle should have similar 
characteristics (such as length, magnitude, etc.). Then we 
make predictions for the current solar cycle 25 and the 
upcoming solar cycle 26. Section 2 describes in detail the data 
and method used in this work. The main results of the long- 
term evolution of solar activity and the prediction of the 
forthcoming solar cycles, as comparisons with the other results, 
are presented in Section 3, and finally, the conclusions of this 
work are summarized in Section 4. 


2. Data and Methodology 


The data used in this work are ASN and the 13 month 
smoothed monthly total sunspot number which are compiled 
and released by the World Data Center (WDC)—Sunspot 
Index and Long-term Solar Observations (SILSO), Royal 
Observatory of Belgium, Brussels. They can be downloaded 
from the website: www.sidc.be/SILSO/datafiles, which 
includes 323 yr of ASN data recorded from 1700 AD to 
present. This is the longest continuous solar activity time series 
available to date, and is one of the most important basic data 
sets for studies of the long-term evolution of solar activity. 

In order to better reflect the long-term trend of the cycle 
maximum of ASN of solar cycles, we first assume a sinusoidal 
empirical function, which can be expressed as 


R= A+ Bsin(22 + Dr) + Ex (1) 


Here x represents the year after 1700 AD, and R means the 
value of ASN corresponding to that year. A, B, C, D and E are five 
unknown parameters that need to be determined by us. The 
assumed empirical function (Equation (1)) is a nonlinear fit 
function, and the parameters are determined in this paper using 
mathematical analysis software. Here, the constant A represents 
level of the background, B is the magnitude of the century cycle, 
C is the period, D signifies the initial phase and E reflects a long- 
term trend of the century cycles. Obviously, the century cycle is 
quasi-periodic. For convenience, we first determine the parameters 
C and D before fitting the function. The wavelet analysis method 
is used to analyze the ASN sequence and to determine whether the 
century cycle ends or not by combining the actual data with the 
consideration of the quasi-periodicity of solar cycles, and finally to 
determine the values of parameters C and D. 

Tan (2011) found that E=0.06, which indicates a slowly 
increasing trend in solar activity on the century timescale. 
Javaraiah (2017) also found that the maximum of solar cycles 
has a linear tendency to increase, and it has also been argued that 
there is no tendency for the maximum of the solar cycle to 
increase since the Maunder Minimum (Svalgaard 2011). There- 
fore, the empirical function assumed in this work also contains the 
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Figure 1. (a) shows the variance plot of the wavelet coefficients after the wavelet transform of the ASN sequence, (b) depicts the contour plot of the real part of the 


wavelet coefficients after the wavelet transform of the ASN sequence. 


aspect of linear variation, and its coefficient is denoted by £. The 
function is fitted to determine whether there is a trend for the 
maximum of the solar cycle to increase (E>0) or 
decrease (E < 0). 

This work assumes that a Schwabe cycle located at the same or 
similar phases in the century cycle will have the same or similar 
cycle characteristics, such as cycle maximum, rise time and cycle 
length. Then, we may adopt the phase similarity method to predict 
the basic features of solar cycles 25 and 26. The specific steps are 
as follows: (1) determine the empirical function corresponding to 
the century cycle, (2) extend the fitting function forward to 
determine the phases of solar cycles 25 and 26 in the century 
cycle, (3) extend the fitting function backward to determine the 


previous cycles that have the same or similar phase with them, (4) 
compare solar cycles 25 and 26 with their corresponding similar 
cycles and determine their characteristic parameters. 


3. Results and Discussion 
3.1. Result 


3.1.1. Determination of the Century Cycle 


Here, we choose the wavelet transform method to analyze 
ASN time series. To ensure the analysis accuracy, we subtract 
the mean from the original ASN time series and then do a 
symmetric extension of the data on both sides to reduce 
boundary effects. 
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Figure 2. The variation of wavelet coefficients in the time domain at 104.0 yr period. 


Figure 1(a) shows the tendency of the variance to change with 
the period. Here, we can find that variance at about a 104.0 yr 
period is the maximum in all periods, which means that the 
104.0 yr period is the most significant periodic component of the 
solar activity. Furthermore, we also find that in terms of peak 
sharpness, it is obvious that the peak of the solar cycle is sharper, 
while the peak of the 104.0 yr period is wider, which indicates that 
the quasi-periodic characteristic of the 104.0 yr cycle is stronger. 
The timescale of the solar cycle obtained by wavelet analysis is 
slightly less than the well-known 11 yr. Due to the quasi-periodic 
characteristic of the solar cycle, we believe that the value is within 
the reasonable range of the timescale of the solar cycle. If we 
assume that the duration values of the 29 complete solar cycles 
observed since 1700 AD conform to a normal distribution and 
perform a Gaussian function fitting, the mean value is around 
10.78, which is close to the value we obtained. The third obvious 
crest is the 56.20 yr period, numerically close to half of the 
century cycle. In order to more intuitively describe the periodicity 
of ASN time series, Figure 1(b) shows the contour plot of the real 
part of the wavelet coefficients of ASN time series. The contour 
plot of the real part of the wavelet coefficients can reflect the 
periodical changes of ASN time series in a different period and its 
distribution in the time domain, and then can be used to judge the 
future trend of ASN in a different period. As can be seen in 
Figure 1(b), there are multiple cycle components in the ASN time 
series during the evolution from 1700 AD to 2022 AD. The 
cyclical components of 10.74 and 104.0yr are consistently 
present, whereas the 56.20 yr period is only partially significant on 
the time axis. This hints at the possibility that the 56.20 yr cycle is 
not a stable periodic feature of solar activity. 

The wavelet coefficients at the 104.0 yr period are plotted in 
Figure 2. We find that there are three full cycles. In addition, 
the absolute amplitude in Figure 2 increases gradually with 


time, which also indicates that there is a gradual trend of 
increasing solar activity since 1700 AD. 

Different works have come to different results on the timescale 
of the century cycle (Gleissberg 1939; Garcia & Mouradian 1998; 
Peristykh & Damon 2003). However, based on our analysis with 
the longest data series so far in this work, we believe that the more 
accurate timescale of the century cycle should be 104.0 yr, which 
is close to the results of 101 yr (Le & Wang 2003), 103 yr 
(Tan 2011) and 100 yr (Feynman & Ruzmaikin 2014), and far 
from the results of 90 yr (Le Mouël et al. 2017). 

After determining the timescale of the century cycle, we 
need more precise confirmation of the initial phase of the 
century cycle on ASN time series to fit the empirical trend 
equation. 

In an earlier study of the century cycle, Tan (2011) suggested 
that solar cycle 24 is located in the valley between the 3rd and 4th 
century cycles. We compare the 13 month smoothed monthly total 
sunspot numbers of the corresponding months in the rising phase 
of solar cycles 24 and 25 one by one. The start of the two solar 
cycles begins with the month in which the minimum of the 
13 month smoothed monthly total sunspot numbers occurs, 2008 
December and 2019 December, respectively. The results are 
shown in Figure 3. We find that the 13 month smoothed monthly 
total sunspot numbers for solar cycle 25, as of 2023 February, are 
overall greater than the corresponding sunspot numbers for solar 
cycle 24. In other words, solar cycle 25 is significantly stronger 
than solar cycle 24. 

Then, how do we determine the end time of the last century 
cycle (G3) and the start time of the forthcoming century cycle 
(G4)? From our above analytical results, we know that solar 
cycle 24 is the weakest cycle in the past half century, and the 
current solar cycle 25 is obviously stronger than solar cycle 24. 
This fact strongly indicates that solar cycle 24 is precisely 
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Figure 3. The red curve presents the smoothed monthly mean sunspot number during the rising phase of solar cycle 25 (from 2019 December to 2023 April), while the 
black curve represents the results during the rising phase of solar cycle 24 (from 2008 December to 2012 April). 


150 
130 
110 
90 
70 
50 
30 


Sunspot number 


0 50 100 150 


200 


Smooth ASN 
—— Fitting fuction 


250 300 350 400 450 


Time/year 


Figure 4. The black curve represents the smoothed ASN series over time, and the red curve is the fitted function of the century cycles. 


located at the lowest point between G3 and G4. Therefore, we 
may conclude that the end of G3 and the beginning of G4 occur 
in the middle of solar cycle 24, that is to say, around 2014 
April. Moreover, solar cycle 25 is already in the rising phase of 
the forthcoming century cycle G4. Going back in time, we can 
see that 1702 AD is the starting year of the first full century 
cycle to appear in the ASN time series. 


3.1.2. Fitting of Hypothesized Empirical Function 


This paper focuses on the long-term trend of solar activity, 
known as the century cycle, and utilizes it to predict the 
parameters of the forthcoming solar cycles. Therefore, in order to 
prevent our study from being affected by outliers in ASN time 
series, we smoothed the ASN time series to better emphasize the 
long-term trend: the maximum of each solar cycle is averaged 
with its minimum at the beginning and with its minimum at the 
end. The resulting averages are used as elements of the smoothed 
ASN time series, and the corresponding times are the closest years 
obtained by indexing real ASN. 


The ASN sequence we used was recorded from 1700 AD, 
but the existence of sunspots was confirmed much earlier 
than 1700 AD. So, in order to minimize the error, this 
paper takes 1700 AD as the starting point, with a total of 
58 sets of data in terms of years, ending at the end of solar 
cycle 24. Its processed image is shown as the black curve 
in Figure 4. 

We assume the empirical function has the form of 
Equation (1), where the parameter C is the length of the 
corresponding century cycle of 104.0 yr, and parameter D, 
which can be calculated from the fact that 1702 is the starting 
year of the century cycle G1, has a value of —0.5385. 

Then, by using the least squares fitting method, we obtain the 
values of parameters A, B, and E in Equation (1), and the 
function can be expressed as 


TX 


R = 78.84 + 23.86 sin 2 
104.0 


— 0.53857) + 0.06028x. 
(2) 
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Table 1 
Summary of Solar Cycle Parameters 
Numbering to Riin Tup Taat Rinax Tae L Rat 
(yr) (yr) (yr) (yr) (yr) 

E 1700 8.300 3 1705 96.70 6 11 0.8330 
D 1711 0.000 6 1717 105.0 6 12 1.000 
C 1723 18.30 4 1727 203.3 6 10 0.6670 
B 1733 8.300 5 1738 185.0 6 11 0.8330 
A 1744 8.300 6 1750 139.0 3 11 1.200 
1 1755 16.00 6 1761 143.2 5 11 1.200 
2 1766 19.00 3 1769 176.8 6 9 0.5000 
3 1715 11.70 3 1778 257.3 6 9 0.5000 
4 1784 17.00 3 1787 220.0 11 14 0.2720 
5 1798 6.800 6 1804 79.20 6 12 1.000 
6 1810 0.000 6 1816 76.30 T 13 0.8570 
7 1823 2.200 7 1830 117.4 3 10 2.333 
8 1833 13.40 4 1837 227.3 6 10 0.6670 
9 1843 18.10 5 1848 208.3 8 13 0.6250 
10 1856 8.200 4 1860 182.2 7 11 0.5710 
11 1867 13.90 3 1870 232.0 8 11 0.3750 
12 1878 5.700 5 1883 106.1 6 11 0.8330 
13 1889 10.40 4 1893 142.0 8 12 0.5000 
14 1901 4.600 4 1905 105.5 8 12 0.5000 
15 1913 2.400 4 1917 173.6 6 10 0.6670 
16 1923 9.700 5 1928 129.7 z 10 1.000 
17 1933 9.200 4 1937 190.6 7 11 0.5710 
18 1944 16.10 3 1947 214.7 7 10 0.4290 
19 1954 6.600 3 1957 269.3 T 10 0.4290 
20 1964 15.00 4 1968 150.0 8 12 0.5000 
21 1976 18.40 3 1979 220.1 7 10 0.4290 
22 1986 14.80 3 1989 211.1 T 10 0.4290 
23 1996 11.60 4 2000 173.9 8 12 0.5000 
24 2008 4.200 6 2014 113.3 3 11 1.200 
25 2019 3.600 


The coefficient of the linear term in the equation is 0.06028, 
which indicates that the intensity of solar activity tends to 
increase linearly on a timescale of centuries. From Figure 4, we 
can also clearly see the presence of the century cycle. 

The fitted function is extended outward by 150 yr, as shown 
by the red curve in Figure 4. Parameters of future solar cycles 
are inferred by the principle of extrapolation of similarity of the 
parameters of solar cycles in the same phase of the century 
cycle. 


3.1.3. Predictions for Solar Cycles 25 and 26 


The method in this paper is a phase-similar prediction based 
on the evolutionary trend of the century cycle, so before the 
prediction, we have to make a collation of the parameters of all 
solar cycles so far. 

The international numbering of solar cycles began in 1755, 
and there have been 24 full solar cycles, so far. In the period 
from 1700 to 1755, ASN data also went through five full solar 
cycles, which for the purposes of analysis we have named 
alphabetically from the latest to the earliest as A, B, C, D and E 
(Tan 2011). 


In addition to the numbering, the characteristic parameters 


used to describe the solar cycle are as follows: 


(a) The year fo of the start of the solar cycle (referred to as the 
start time); 

(b) Minimum of ASN at the start of the solar cycle, Rmin; 

(c) The length of the rising phase of the solar cycle, Typ 
(referred to as the rise time); 

(d) The year of the maximum of the solar cycle, tmax (referred 
to as the peak time); 

(e) Maximum of ASN of the solar cycle, Rmax; 

(f) The length of the decay phase of the solar cycle, Ty. 
(referred to as the fall time); 

(g) The period of the solar cycle, L (referred to as the period); 

(h) Asymmetric parameter, Ra, the ratio between the rise 
time and the fall time, which represents whether the shape 
of the solar cycle is left-right symmetrical or not. 


The parametric characteristics of Schwabe cycles are listed 


in Table 1. 


Table 1 lists the parametric characteristics of 29 solar 


Schwabe cycles since 1700 AD. First of all, the period L takes 
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Figure 5. The black curve represents the trend of the observed ASN with time, and A-E and 1-26 are numbers of the Schwabe solar cycles. The red curve represents 
the fitted curve of the century cycles, where G1, G2, G3, and G4 number the different century cycles, respectively. The blue curve is the predicted profiles of the 


forthcoming solar cycles 25 and 26. 


values between 9 and 14 yr, and 82.76% of the period L of all 
recorded solar cycles is between 10 and 12 yr. This confirms 
the quasi-cyclical nature of the solar cycles. Second, there is a 
big difference among the maxima of the different solar cycles. 
So far, the smallest observed maximum occurred during solar 
cycle 6 with a maximum ASN of 76.30. The highest maximum 
occurred during solar cycle 19 with a maximum ASN of 269.3. 
In addition, the following conditions are observed in the 
asymmetric parameter R}: 75.86% of Rx of solar cycles are less 
than 1, 13.79% of Rw of solar cycles are greater than 1 and 
10.34% of Ra of solar cycles are equal to 1. This indicates that 
most solar cycles are rapidly rising and slowly falling, 
belonging to the left-biased asymmetric type. 

After clarifying the main properties of the solar cycles, we 
try to make a prediction for the forthcoming solar cycles 25 
and 26. 

First, let us determine the phase of the upcoming solar cycles 
25 and 26 in the century cycle and other Schwabe cycles with 
similar phases. Figure 5 plots the temporal profiles of 29 
Schwabe cycles numbered with 1—24 and A-E and overlaid on 
the century cycles numbered G1-—G4 since 1700, respectively. 
It shows that solar cycle 24 is located at the valley between the 
century cycles G3 and G4, while solar cycles 25 and 26 are the 
first and second cycle in the ascending phase of the century 
cycle G4. Considering the solar cycles located in the early part 
of the century cycle G1 (A, B, C, D, and E), for which there is a 
large amount of missing data (Clette et al. 2015; Clette & 
Lefévre 2016), we do not consider solar cycles located in the 
early part of the century cycle G1 when searching for the 
corresponding past solar cycles. 

We find that the phase of solar cycle 25 is very similar to that 
of solar cycle 6 in century cycle G2 and solar cycle 15 in 
century cycle G3, while the phase of solar cycle 26 is very 
similar to that of solar cycle 7 in century cycle G2 and solar 
cycle 16 in century cycle G3. Naturally, we may apply the 


averaged parameters’ values of solar cycles 6 and 15 to predict 
the main characteristics of solar cycle 25, and apply the 
averaged parameters’ values of solar cycles 7 and 16 to predict 
the main characteristics of solar cycle 26. At the same time, we 
have to consider the effect of the last linear increasing term in 
Equation (2). 

The details of our prediction are as follows: 

About solar cycle 25: 


1. Rmax Of solar cycle 6 is 76.30, and Rmax Of solar cycle 15 
is 173.6. Based on the principle of similarity, we can 
speculate that Rmax of solar cycle 25 should be between 
76.3 and 173.6. In addition, since the corresponding 
value of solar cycle 25 in the phase of the century cycle is 
higher than that of solar cycle 24 in the trough, Rmax of 
solar cycle 24 is 113.3, and there is an overall 
enhancement trend over the centuries, so we finally 
predict that Rmax of solar cycle 25 will be between 113.3 
and 180.1, that is to say: 146.7 + 33.40; 

2. L for solar cycle 6 has a period of 13 yr and solar cycle 15 
has a period of 10 yr. In previous analyses of the overall 
characterization of solar cycle parameters, long cycles of 
13 yr accounted for only 6.9% of all cycles, including 
solar cycle 6 and solar cycle 9. Solar cycle 6, which 
began in 1810, is in the historically famous Dalton 
Minimum (1790 AD-1830 AD). Based on Table 1 and 
the century cycle phase, the intensity of solar cycle 25 is 
not predicted to be lower than that of solar cycle 24. 
Therefore, we tentatively believe that solar cycle 25 will 
not enter a new minimum similar to the Maunder 
minimum and the Dalton minimum. Regarding the 
relationship between the period of solar cycle and the 
grand minimum, it was pointed out in a study (Karak & 
Choudhuri 2013) that the period may be prolonged in 
1-2 solar cycles before entering the grand minimum. 
Combining these analyses, we believe that the probability 
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of a 13 yr length for solar cycle 25 is very low, and 
therefore predict that its period L will be about 11 + 1 yr. 

3. Rat of solar cycle 6 and R, of solar cycle 15 are both less 
than 1, so we predict that R,, of solar cycle 25 is also less 
than 1, i.e., the rise time is less than the fall time. 

4. Typ, with the predictions of the period and asymmetric 
parameter, the rise time predicted in this paper is 4—5 yr. 
Since the actual start of solar cycle 25 is 2019 December, 
we predict that Rmax Of solar cycle 25 will occur in 2024. 


About solar cycle 26: 


1. Rmax Of solar cycle 7 is 117.4, Rmax Of solar cycle 16 is 
129.7. we apply the averaged value of solar cycle 7 and 
16 to predict solar cycle 26. Considering the increasing 
trend of solar cycles (the last term of Equation (2)), we 
finally predict that Rmax of solar cycle 26 should be 
133.0 + 3.200; 

2. L, the period of both solar cycle 7 and solar cycle 16 is 
10 yr, so the period of solar cycle 26 is also predicted to 
be 10 yr; 

3. Rat of solar cycle 7 is greater than 1 and Ry of solar cycle 
16 is equal to 1. Therefore, Ra of solar cycle 26 is 
predicted to be greater than or equal to 1, i.e., the rise 
time will not be shorter than the fall time; 

4. Typ, the rise time of solar cycle 7 is 7 yr and the rise time 
of solar cycle 16 is 5 yr, because R(= 2.333) of solar 
cycle 7 is the highest of all solar cycles. Considering the 
period L and asymmetric parameter R, this paper 
predicts a rise time of 5—6 yr for solar cycle 26. In the 
previous section, our prediction for the period L of solar 
cycle 25 is 11 + 1 yr, and in predicting solar cycle 26, we 
take the average of the predicted values for L of solar 
cycle 25. In other words, we consider the end of solar 
cycle 25 to be 2030, and the beginning of solar cycle 26 
to be 2030. Therefore, we predict that Rmax of solar cycle 
26 will occur between 2035 and 2036; 

5. For Rmin, according to the fitted empirical function, there 
is an overall tendency for solar activity to increase, so we 
predict that Rmin of solar cycle 26 will be greater than or 
equal to 3.6 of solar cycle 25. 


We plotted the prediction curves for solar cycles 25 and 26 
as shown by the blue dotted curve in Figure 5. The values are 
obtained by averaging the previous two solar cycles plus the 
appropriate solar activity trend. It should be noted that our 
predictions above do not take into account the influence of the 
56.02 yr period. Although the 56.02 yr period is the third 
strongest in Figure |(a), its presence is not stable (Figure 1(b)), 
and we have reason to suspect that solar cycle 25 is not affected 
by this period. However, if we consider the influence of the 
56.02 yr period, we can also see from Figure 5 that the first 
three century cycles all fluctuate at their peaks. The troughs are 
located in solar cycle A and there is one in the century cycle 
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G1, solar cycle 10 in the century cycle G2, and solar cycle 20 
in the century cycle G3. However, solar cycles 25 and 26 are 
located in the ascending phase of the century cycle, so this does 
not affect our predictions. 

The good fit of the actual ASN sequence to the fitted century 
cycle curve in Figure 5 also further validates the correctness of 
the timescale of the century cycle and the century cycle G3 
ending times that we identified. 


3.1.4. Comparison with Other Predictions 


In addition to the prediction of solar cycle 25, this paper also 
provides a new prediction for solar cycle 26, which we will 
compare with the results of the existing predictions of solar 
cycles 25 and 26 respectively in the following. 

Up to now, many results for solar cycle 25 predictions have 
been published. Nandy (2021) has done a detailed compilation 
of most of the solar cycle 25 predictions prior to the start of 
solar cycle 25, and since the predictions of solar cycle 25 in this 
paper are already in the ascending phase of solar cycle 25, we 
have chosen more predictions made after the beginning of solar 
cycle 25 (after 2019). We summarize the results of the recent 
predictions, sorted by chronological order of the predictions in 
the different categories of methods (as shown in Table 2). It 
should be noted in Table 2 that there are also differences in the 
accuracy of the results obtained by different prediction 
methods, which is also a manifestation of the characteristics 
of different methods, so we do not make a uniform requirement 
for accuracy here. Second, according to the classification of 
prediction methods in the previous section, our prediction 
method belongs to the extrapolation method, which is listed 
separately in the table. 

Comparing the Rmax of all predictions in Table 2, we can 
see that: 


(a) Guo et al. (2021) in the dynamo model category, Diego & 
Laurenza (2021), Bisoi et al. (2020) in the precursor 
method, McIntosh et al. (2020), Sarp et al. (2018) in the 
extrapolation method, and Su et al. (2023), Prasad et al. 
(2022) in the interdisciplinary category are consistent 
with our results that solar cycle 25 will have a larger Rmax 
than solar cycle 24. Among these results, Diego & 
Laurenza (2021), who utilize the relationship between the 
recurrence index of geomagnetic activity in the falling 
phase of the previous cycle and sunspot number in the 
next cycle for prediction, reach the same conclusion as 
ours that the century cycle ends at solar cycle 24; 

(b) Nagovitsyn & Ivanov (2023), Lu et al. (2022), Xiong 
et al. (2021), Kumar et al. (2021) in the precursor method, 
and from their predicting results for Rmax, We also believe 
that their conclusion is that solar cycle 25 is stronger than 
solar cycle 24; 

(c) In addition, as in Javaraiah (2023) of the precursor 
method, Javaraiah (2023) of the extrapolation method and 
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Table 2 


Comparison of Different Predictions for Solar Cycle 25 
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Forecast Information 


The Prediction Result of Solar Cycle 25 


Category Serial Number Time of Publication Author Rigg tae 
The Century Cycle a) 2023 This paper 146.7 + 33.38 2024 
Dynamo Model Category (2) 2021 Guo et al. 126 e 
(3) 2019 Labonville et al. sotit 2025.37 083 
Precursor Method (4) 2023 Nagovitsyn & Ivanov 149 + 28 2023.5—2024.5 
(5) 2023 Javaraiah (polar field) 1257 te 
(6) 2022 Brajsa et al. 121 + 33 tee 
(7) 2022 Lu et al. 145.3 2024 Oct 
(8) 2021 Burud et al. (spotless days) 99.13 + 14.97 2024 Feb-—2024 Mar 
(9) 2021 Burud et al. (aa index) 104.23 + 17.35 2024 Feb-2024 Mar 
(10) 2021 Diego & Laurenza 205 + 29 around mid-2023 
(11) 2021 Xiong et al. 140.2 2024 Mar 
(12) 2021 Chowdhury et al. 100.21 + 15.06 2025 Apr + 6.5 months 
(13) 2021 Kumar et al. 126+3 te 
(14) 2020 Bisoi 133 + 11 
Extrapolation Method (15) 2023 Javaraiah (sunspot group area) 125+ 11 
(16) 2021 Kakad & Kakad 103.3 + 15 
(17) 2020 McIntosh et al. 229 + 76 te 
(18) 2018 Sarp et al. 154 + 12 2023.2 + 1.1 
Other Category (19) 2021 Courtillot et al. 97.6 + 7.8 2026.2 + 1 
(20) 2019 NOAA/NASA 115+ 10 2025 Jul + 8 months 
Interdisciplinary Category (21) 2023 Su et al. 133.9 + 7.2 2024 Feb 
(22) 2022 Bizzarri et al. (the TT model) 104+7 2024 Jul + 7 months 
(23) 2022 Bizzarri et al. (the CV model) 110+9 2024 Apr + 12 months 
(24) 2022 Prasad et al. 171.9+3.4 2023 Aug + 2 months 
(25) 2020 Attia et al. 80 + 12 2026 


the result (20) which is published by the NOAA/NASA 
co-chaired international panel to forecast Solar Cycle 25 
(https: //www.swpce.noaa.gov /news /solar-cycle-25-forecast- 
update), these predictions of Rmax are similar to Rmax Of solar 
cycle 24; 

(d) Labonville et al. (2019) in the dynamo model category, 
Brajša et al. (2022), Burud et al. (2021), Chowdhury et al. 
(2021) in the precursor method, Kakad & Kakad (2021) 
in the extrapolation method, Courtillot et al. (2021) in the 
other category and Bizzarri et al. (2022), Attia et al. 
(2020) in the interdisciplinary category, lead directly to 
the conclusion that solar cycle 25 will have a smaller 
Rmax than solar cycle 24, or that cycle 25 lies in the 
trough between century cycles, or that solar cycle 25 will 
enter a new grand minimum. This is inconsistent with our 
prediction using the century cycle. 


A comparison with the prediction results for cycle 26 is 
shown in Table 3. 

By comparison, we find that the prediction results of Rmax of 
solar cycle 26 are quite different in different studies. Kalkan 
et al. (2023) using nonlinear autoregressive exogenous neural 


networks are considered to be similar to or weaker than solar 
cycle 24. Solar cycles 25 and 26 are believed to be at the 
minimum of the century cycle, and solar activity will enter a 
new grand minimum. Liu et al. (2023) using the long short- 
term memory method show that solar cycles 25 and 26 are 
similar, which is different from our result that solar cycle 26 is 
slightly smaller than solar cycle 25. Becheker et al. (2023) 
using the function fitting method considers solar cycle 26 to be 
numerically weaker than solar cycle 24. Current projections 
generally agree that solar cycle 26 will peak between 2035 
and 2036. 


3.2. Discussion 


The ASN sequence is the longest time series data that 
humans have continuously observed of the same celestial body 
to date, but compared with the history of solar activities its data 
are still too few and too short. Therefore, it is practically 
impossible to tell whether the quasi-periodic properties of the 
solar cycle and the century cycle are due to the incomplete 
observed cycles showing quasi-periodic properties because of 
the short time being observed, or whether they are the 
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Table 3 
Comparison of Different Predictions for Solar Cycle 26 


The Prediction Result of 


Forecast Information Solar Cycle 26 


Serial Time of 
Number Publication Author Rmax fax 
(1) 2023 This paper 133.0 + 3.200 2035-2036 
(26) 2023 Kalkan et al. 113.25 2036 Oct 
(27) 2023 Liu et al. 135 2035 Jan 
(28) 2023 Becheker 96 +28 

et al. 


properties of the solar cycle and the century cycle themselves. 
To solve this problem, it is necessary to study the forward 
extension of the sunspot number sequence. The impact of this 
problem in this paper is that in terms of prediction error, our 
predictions depend on the peak of the solar cycle at a similar 
phase of the century cycle to the one being predicted, but there 
are only two similar-phase solar cycles available for prediction, 
and so this shows a more different range of error in the 
predictions for solar cycles 25 and 26. Our method analyzes 
long-term trends (centennial scales) in solar activity and is 
therefore dependent on the length of the data. This is in contrast 
to the vast majority of prediction methods, such as the 
precursor method considering polar magnetic field, which 
utilizes only the polar magnetic field of the previous cycle on a 
scale of about a decade. 

The error in the prediction results obtained by different 
prediction methods varies considerably. For example, the 
method (McIntosh et al. 2020) which has the dependence of 
prediction results on terminators gives a large error range when 
there are undetermined terminators. However, the prediction by 
using comprehensive precursors and multiple regression 
techniques can significantly improve the prediction accuracy, 
adaptability, and stability, and the regression coefficient can 
reach 0.95 (Xiong et al. 2021; Lu et al. 2022). As for the 
prediction results of the interdisciplinary (specifically referring 
to machine learning) category, the goodness of the results lies 
in the selection and construction of the model, so the error 
ranges of the prediction results are generally small after the 
model is determined. However, the quasi-periodicity and 
complexity of the solar activity make the predictions challen- 
ging, and too precise predictions might have a risk of 
overfitting. 

In fact, our prediction also includes the effect of the 56.20 yr 
period component in the Rmax of the solar cycles. This is 
because we assume that cycles located in similar phases on the 
background of the century cycle will have similar parameter 
characteristics, while the 56.20 yr period component is only 
implied in the relatively small Schwabe cycle of the peak year 
of the century cycle. Figure 5 shows that cycle A in G1, cycle 
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10 in G2 and cycle 20 in G3 are located near the peak year of 
the century cycles, but they are similarly relatively weak. Based 
on their statistical trend, we may also predict the future cycle 
located near the peak year of G4. Obviously, this component 
has little impact on the prediction of solar cycles 25 and 26 
which are located in the early ascending phase of the century 
cycle. 

With regard to the comparison of results, according to 3.1.4 
we can easily realize that the published predictions are all very 
significantly different from each other, presenting a chaotic 
state. We believe that this is a clear manifestation of the 
complex and quasi-periodic nature of the solar cycle. Although 
there is a huge difference between the results, a comparison of 
the results obtained from the different methods of prediction 
can also give us a deeper understanding of the solar cycle. We 
also find that, of all the predictions, the later the prediction, the 
greater the probability of the result of the work that believes 
that solar cycle 25 will be stronger than solar cycle 24, which is 
also consistent with our predictions. One of the reasons for this 
is that the later the prediction is made, the more pronounced the 
evolutionary trend exhibited by the relative sunspot number in 
solar cycle 25, but it also illustrates that none of the current 
types of prediction methods have a clear advantage in 
predicting a long time in advance, with the most constrained 
being the precursor method. Unlike these methods, our 
prediction method is based on the idea that solar cycles in 
similar phases of the century cycle have similar cycle 
parameters. Therefore, after determining the existence of the 
century cycle, as well as the shape (fitting empirical function) 
and timescale of the century cycles (parameter C), our 
prediction is not basically limited by the early or late prediction 
time. The closer the time of our prediction is to the time at fmax» 
the closer our prediction will be to the true values, without 
affecting our prediction of the long-term trend of the solar 
cycle. Our method has a requirement for data length, but there 
are two sides to everything. Regarding the earlier prediction, 
our work is superior. 

Among the many results of predictions for solar cycle 25, 
many thought solar cycle 25 would be smaller than solar cycle 
24, contrary to our prediction, and they thought that there 
would be a new grand minimum in solar cycle 25 similar to the 
Maunder minimum and the Dalton minimum. However, 
looking at the current data for solar cycle 25, the probability 
of it being weaker than solar cycle 24 is dramatically reduced. 
In other words, the trough between century cycles G3 and G4 
does not have a very pronounced long-lasting grand minimum, 
and we can probably assume that the correlation between 
century cycles and grand minimum is weakened. But how 
exactly this will actually play out, and whether a new grand 
minimum will emerge, still depends on how the sunspot 
number data develop. 
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4. Conclusion 


This work confirms the existence of solar century cycles, 
determines their main parametric characteristics, and obtains an 
empirical function (Equation (2)). 


1. The century cycles do exist and their period L is 104.0 yr. 

2. Solar activity has a gradually enhancing trend on the 
timescale of century cycles. 

3. The solar Schwabe cycle 24 is located in the valley 
between century cycles G3 and G4. After that, cycles 25 
and 26 are in the ascending phase of century cycle G4, 
therefore, a new grand minimum will not occur. Or, that 
is to say, the grand minimum just occurred around 
cycle 24. 


The existence of the solar century cycle indicates that there 
must be a global physical mechanism on the Sun that is larger 
in scale than the solar Schwabe cycle, influencing the patterns 
of solar activity. But so far, we are not yet understanding this 
global physical mechanism. The century cycle pattern can help 
us predict the basic characteristics of forthcoming Schwabe 
cycles. Here, based on the empirical functions of the solar 
century cycle, and the assumption of similarity extrapolation, 
we predict the main parametric characteristics of solar cycles 
25 and 26. The prediction results are as follows: 


1. Solar cycle 25 will reach its maximum in 2024, with a 
peak ASN of 146.7 + 33.40, stronger than cycle 24, and 
its period is about 11 + 1 yr; 

2. Solar cycle 26 will start from 2030, and reach its 
maximum between 2035 and 2036 with a peak ASN of 
133.0 + 3.200, and its period will be about 10 yr. 


This paper ignores some sunspot anomaly data when 
examining the century cycle of solar activity. This may result 
in our predictions of the long-term trend not being wrong, but 
there are some uncertainties in the details. It is possible that 
future studies will be able to interpret the sunspot anomaly data 
and combine multiple solar cycle features for future cycle 
predictions to improve accuracy. In addition, in studying the 
long-term evolutionary characteristics of solar activity, long- 
term data on isotopes (especially carbon-14 isotopes) on Earth 
could be analyzed in addition to observations of the relative 
sunspot number. This has the potential to extend the time series 
to thousands of years in length and give a more accurate picture 
of the solar activity’s long-term evolutionary pattern. 
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